Static impurities in a supersolid of interacting hard-core bosons on a triangular lattice 
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We study the effect of impurities in a supersolid phase in comparison to the behavior in the solid 
and superfluid phases. A supersolid phase has been established for interacting hardcore bosons 
on a triangular lattice which may be realizable by ultracold atomic gases. Static vacancies are 
considered in this model which always lower the magnitude of the order parameter in the solid 
or superfluid phases. In the supersolid phase, however, the impurities directly affect both order 
parameters simultaneously and thereby reveal an interesting interplay between them. In particular 
the solid order may be enhanced at the cost of a strong reduction of the superfluidity, which shows 
that the two order parameters cannot be in a simple superposition. We also observe an unusual 
impurity pinning effect in the solid ordered phase, which results in two distinct states separated by 
a first-order transition. 
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A bosonic supersolid phase is characterized by the co- 
existence of two seemingly contradictory order param- 
eters, a solid crystalline order and a superfluid den- 
sity. This reflects the spontaneous breaking of two in- 
dependent symmetries, namely translation and a U(l) 
gauge rotation, which are also known as diagonal and off- 
diagonal order, respectively. The simultaneous breaking 
of two independent symmetries in the supersolid phase is 
counter-intuitive and unusual, because normally a spon- 
taneously broken order locks the system into a single 
phase. Only when the remaining fluctuations are large 
enough, two independent order parameters may exist in 
one phase, e.g. due to frustration. Having been predicted 
forty years agoji supersolids recently received renewed 
interest after a possible observation in 4 Ke£ The pres- 
ence of mobile 3 He impurities appears to be important 
in those systems, which are predicted to raise T c but re- 
duce the superfluid density^ While the experiments are 
still controversially discussed,— there is now very strong 
numerical evidence that a supersolid phase is realized for 
interacting hardcore bosons on a triangular lattices- 
Such a model can potentially be realized by ultracold 
atoms in optical trapsJ^ Tunable superexchange models 
have already been experimentally created in this rapidly 
emerging field^ and interacting hard core bosons^ have 
also been discussed. Most recently also triangular lattices 
have been possible^ so it appears likely that supersolid- 
ity will soon be a central topic for hard core boson models 
in the ultra-cold gases community. 

While the coexistence of the corresponding two order 
parameters is well established numerically in the hard 
core boson systems^ the microscopic interplay between 
them is still unclear. We now study the impurity effects 
on both order parameters simultaneously in the super- 
solid phase, in order to clarify if the two order parameters 
are in a simple superposition or how they may interact 
locally. The use of substitutional impurities in strongly 
correlated systems has become a standard tool for un- 
derstanding the underlying quantum phases In par- 



ticular, it is possible to study local expectations values 
around defect o 20 i 21 for an analysis of the elementary ex- 
citations and direct comparison with theoretical models. 
In the supersolid phase we are now able to consider the 
effect of static impurities on two coexisting order param- 
eters simultaneously. 

The model we will consider in this paper is the spin-1/2 
model on a triangular lattice 

H = -tJ2(S+Sr + h.c.) + VJ2S;S] -BJ2^, (1) 

<y> (ij) » 

with antifcrromagnctic exchange V in the z-direction cor- 
responding to nearest neighbor repulsion and ferromag- 
netic exchange t in the x-y-directions, corresponding to 
the kinetic energy and B = fx — 3V in terms of the chem- 
ical potential of the equivalent hard-core boson problem. 
The simplest impurities are given by lattice vacancies in 
model (JTJliS 2 . 

The two order parameters in the supersolid arc 
given by the structure factor S(q = [4/37T, 0]) = 
(I J2 r Sr e%q ' r \ 2 )/N 2 for the solid order, and by the su- 
perfluid density p s , which is typically measured using 
the winding number W in quantum Monte Carlo (QMC) 
simulations p s = (W 2 ) /4ffl£ ~ 11 i 23 We use a modified per- 
turbation theory and the directed-loop stochastic series 
expansion QMC algorithm 2 ^ with finite size scaling up 
to N = 324 sites at a temperature of T = 0.02V. In 
order to avoid trapping in one of the degenerate states it 
is also essential to implement parallel tempering in the 
parameter space i 25 ' 26 

For reference we first examine a single vacancy in the 
solid phase, which already shows interesting effects. As 
shown in the inset of Fig. [T] this phase occurs for small 
xy-coupling t. The solid order is characterized by a 2/3 
filled state for positive fields with exactly two spins on 
each triangle pointing up. For negative fields there is 
an equivalent 1/3 filled ordered phase due to the spin 
flip symmetry around B = 0. The vacancy in the XXZ 
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FIG. 1: (color online). The average magnetic moment of the 
sublattice which contains the vacancy (solid) and the other 
two sublattices (dashed). Inset: The phase diagram with the 
impurity phases and the plot trajectory we use in Figs. [TJ [3] 
and 03 along t/V = 0.08 and then along B/V = 0.2. 



model U) docs not break this symmetry,— so it is suffi- 
cient to consider only positive fields B > in the phase 
diagram in Fig. [T] The choice of the spin-down sublattice 
(pointing against the field) gives a three-fold degeneracy, 
which is however lifted by the vacancy. In particular, for 
< B < 1.5V and t = the vacancy site must belong 
to one of the spin-up sublattices, as can be seen by sim- 
ple energetic considerations. Therefore the order in the 
entire system is pinned by a single defect and only a two- 
fold degeneracy remains. The average occupation (S z ) on 
the different sublattices in Fig. [T] shows that this pinning 
also continues throughout the supersolid phase. How- 
ever, the spin density of the other two sublattices sur- 
prisingly point against the field in the supersolid phase. 

For larger fields B > 1.5V there is a transition to a 
different state, where the order is now pinned on the op- 
posite sublattice with no remaining degeneracy. There- 
fore, a single impurity can in fact induce a transition 
between two distinct ordered states of the entire sys- 
tem. The transition line also depends on the xy-coupling 
t as can be seen by perturbation theory in the "hop- 
ping" terms Hij = —t(S^Sj + h.c). Unfortunately, 
the usual perturbative correction to the wavefunction 
IV) « |0> + £<«> lyXijlfly-IOV^o - Eij) diverges with 
the number of hopping terms, i.e. the system size N. 
Here |0) is the ordered state and has opposite spins 
exchanged on the bond i,j relative to |0). Of course the 
number of lattice sites N must be irrelevant in the or- 
dered phase, so the trick is to modify the perturbation 
correction to include only those hopping terms which ac- 
tually affect a local expectation valued For example, to 
calculate the energy correction at one bond 8E+j only the 
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FIG. 2: (color online), (a) {SfSj §1} with i,j,k on one tri- 
angle, (b) (SfSj) and (c) {S/~S~~}, where i,j are neighboring 
sites on sublattices not occupied by the impurity as a func- 
tion of distance in the solid phase (t/V = 0.08, B/V = 0.65, 
N — 144). The modified perturbation theory (solid) agrees 
very well with the QMC results (dashed). Inset: Distribution 
of S7) and (SiSjSk) on the bonds and triangles, respec- 
tively. Red signals an increase and blue a decrease relative to 
the gray bulk values. 



corresponding hopping term is considered 

Wa) « io> + imm^/iEo - E i:j ) (2) 

and we simply get 

|(ij|i^|0}| 2 /(£ - Ey) = t 2 /(E - E i:i ). (3) 



For the case of a vacancy, the excitation energy Eij de- 
pends on the location of the bond and it also depends 
on which sublattice is pointing down in the ordered state 
0). After summing over all contributions, we find that 
the energy difference between the two possible pinned or- 
dered states is given by AE = B - ^ - ^ + 0(t 3 ). 
Therefore, the impurity driven transition line runs along 



(4) 



as shown in the inset of Fig. [T] which also agrees with 
our numerical QMC results. 

Using the modified perturbation theory with a re- 
stricted sum in Eq. ^ it is also possible to analytically 
calculate local expectation values, e.g. when calculating 
(S?) all hopping terms connecting to the site i are in- 
cluded. The results give a good indication about the lo- 
cal order around the vacancies. In particular, (SfSj) and 

(SfSjSf.) with i,j,k on neighboring sites are good indi- 
cators of the local solid order, which are reduced around 
the vacancy. On the other hand the quantum fluctuations 
(S^~S~) on bonds are enhanced as shown in Fig. O In- 
terestingly, the order reduction is not correlated in space 
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FIG. 3: (color online). Scaled impurity corrections to the 
superfluid density NAp s and the structure factor NAS(q). 
Inset: Order parameters in a pure system. 



with the fluctuation enhancement and the effect is also 
not always strongest directly at the vacancy. The inset in 
Fig.[Ushows the distribution of the local expectation val- 
ues around the impurity on the lattice directly. The ex- 
pectation values (S^S~) between two non-impurity sub- 
lattices is different from those bonds involving the impu- 
rity sublattice even very far from the impurity. This is 
a secondary effect from the pinned order and should not 
be mistaken for an independent bond order. 

It is not surprising that an impurity generally reduces 
the order parameter locally. The main question for the 
supersolid phase is now if the vacancy reduces both order 
parameters as may be expected for a simple superposi- 
tion of the two effects or if an interesting interplay can be 
observed. The answer to this question is summarized in 
Fig. |31 where we plotted the impurity contributions of the 
two relevant order parameters in the system as we cross 
the phase boundaries along the trajectory in the inset of 
Fig. [TJ The dominant parameter in the phases with one 
single order is always reduced, while the other param- 
eter remains unchanged close to zero. However, in the 
supersolid phase only the superfluid density is strongly 
reduced, while the solid order is in fact enhanced. It is far 
from obvious why the vacancy should enhance the solid 
order in this case, contrary to what we observed in the 
solid phase. The only explanation of the observed be- 
havior is that the vacancy reduces the superfluid density 
locally so strongly that the solid order is revived, which is 
evidence for a microscopic competition between the two 
order parameters. This result clearly shows that the two 
order parameters are not in an independent superposi- 
tion. It is noteworthy that the competition and the total 
change of the order parameters is strongest close to the 
second order phase transitions to the superfluid phase. 

The local expectation values in Fig. @] also demon- 
strate the competition of order parameters locally. The 
local solid order (SfSj-) and close to the im- 
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FIG. 4: (color online). Local expectation values as a function 
of geometrical distance analogous to Fig. [2] in the supersolid 
phase (t/V = 0.08, B/V = 0.2, N = 144). Closest to the 
impurity the values almost recover the bulk values in the solid 
phase in Fig. [2] Inset: Analogous to inset O but showing the 
opposite behavior on a relative color scale. 



purity is now enhanced while the kinetic energy (Sj~S~) 
is strongly reduced. This is in strong contrast to the ob- 
servations in the solid phase in Fig. [5] and the relative 
changes are also much more dramatic and correlated in 
space, which again demonstrates the direct interplay be- 
tween both order parameters. 

We finally turn to the interesting case of several impu- 
rities in the system. A second vacancy on the same sub- 
lattice is rather trivial, corresponding to the same pinned 
order, i.e. constructive interference of the induced mag- 
netization density. Remarkably also a second impurity 
on an opposite sublattice has the same effect since the 
second impurity simply lifts the remaining two-fold de- 
generacy exactly in such a way that both impurities are 
located on spin-up sublattices. The order is now com- 
pletely pinned, but all observed effects are approximately 
additive, i.e. the impurity contributions to the order pa- 
rameters simply double in the entire parameter space 
in Fig. [SJ Also the observed phase transition between 
pinned states in Fig.[T]remains unchanged. This situation 
is in sharp contrast to two impurities on different sublat- 
tices in an unfrustrated square lattice^ which must show 
destructive interference of the alternating magnetization. 

The most interesting case is given by three impurities 
on the three different sublattices. The threefold degen- 
eracy of the lattice is then approximately restored again 
without pinning and the impurities must interfere de- 
structively. This results in a surprisingly strong reduc- 
tion of the solid order S(q) which now also carries over 
into the supersolid phase as shown in Fig. [5j while also 
the reduction of p s remains strong. Obviously the effects 
arc not simply additive in this case and indicate some 
interesting impurity-impurity interactions. However, in 
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FIG. 5: The impurity contributions to the structure factor 
and the superfluid density for two and three impurities on 
different sublattices in a system of 144 sites. 



generic systems impurities break the symmetry between 
the three sublattices, so that the observed pinning and 
order parameter competition described above is the more 
general scenario. 

In summary, we have used a modified perturbation the- 
ory and QMC simulations to analyze impurity effects in 
a supersolid in comparison to other phases with single 
order as realized by the model in Eq. ([1}. 

In the solid phase a non-trivial pinning of the entire 
order by a single defect has been observed. Therefore, 
impurities create a first order transition line between two 



different pinned states given by Eq. (j4|) , which is not seen 
in the pure system. 

In the supersolid phase the solid order is surprisingly 
enhanced by an impurity, which coincides with a strong 
reduction of the superfluid order. This is evidence for 
an interesting microscopic competition between the two 
order parameters, which certainly cannot be in a simple 
superposition. 

For two impurities a simple addition of the observed 
effects can be seen, while for three impurities on different 
sublattices a strong destructive interference changes the 
physics completely. For a more complete understanding 
of the impurity-impurity interactions on three different 
sublattices more research is needed. 

In all impurity configurations a very strong reduction 
of the superfluid density p s occurs close to the second or- 
der supersolid-superfluid transition. It is therefore likely 
that the superfluid order can be destroyed with a critical 
density of impurities, while the solid order may be en- 
hanced. The extreme limit of this effect corresponds to 
the removal of one sublattice, i.e. the honeycomb lattice, 
which indeed results in an extended solid phased! 
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